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We investigate simplified dark matter models where the dark matter candidate is a Dirac 
fermion charged only under a new gauge symmetry. In this context one can understand 
dynamically the stability of the dark matter candidate and the annihilation through the new 
gauge boson is not velocity suppressed. We present the simplest Dirac dark matter model 
charged under the local B — L gauge symmetry. We discuss in great detail the theoretical 
predictions for the annihilation into two photons, into the Standard Model Higgs and a 
photon, and into the Z gauge boson and a photon. Our analytical results can be used for 
any Dirac dark matter model charged under an Abelian gauge symmetry. The numerical 
results are shown in the B — L dark matter model. We discuss the correlation between the 
constraints on the model from collider searches and dark matter experiments. 
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I. INTRODUCTION 

The existence of dark matter (DM) in the Universe has motivated the particle physics community 
to investigate extensions of the Standard Model (SM) of particle physics, with countless dark 
matter candidates on the market which need to be probed against experiments. There are basically 
three ways to search for these candidates: one looks for possible signals in dark matter direct 
detection experiments, for ganrma-ray lines and other signals from dark matter annihilation in 
indirect detection experiments, and for large missing energy signatures together with a mono-jet or 
a mono-photon at colliders. See Refs. [1-4] for a detailed discussion of these possibilities. 

Most of the effort so far has been focused on studies of complete models and their dark matter 
candidates, as for example the neutralino in the Minimal Supersymmetric Standard Model. Also 
effective field theory (EFT) has been used to constrain the scale of the new physics, particularly 
using the data from the LHC run one. However, given the large center-of-mass energy of the LHC 
in its second run, it is obvious that the EFT approach is prone to fail in large fractions of parameter 
space. Therefore, it is sensible to consider simplified dark matter models which can capture the 
main features of the dark matter sector with a limited number of parameters. For a discussion of 
simplified models see Refs. [5, 6]. 

A simplified model should contain a mediator and a dark matter candidate, and it should not 
violate generically any low energy observables. These criteria are necessary but not sufficient, since 
a naive model can result in misleading statements. We argue that the simplified model must be 
self consistent, i.e., not violate gauge invariance and be anomaly free. Only then is it possible to 
perform a full study and compute for example higher order processes leading to gamma-ray lines 
relevant for indirect dark matter searches. 

In this article we investigate in detail simplified models for Dirac dark matter, which have the 
following features: 

• The dark matter is charged under a local gauge symmetry. This symmetry can be sponta¬ 
neously broken at the low scale and a remnant discrete symmetry guarantees the dark matter 
stability. In the simplest case one has a local U(l)' symmetry broken to a Z 2 discrete sym¬ 
metry. I 11 the case of an unbroken gauge symmetry one can use the Stueckelberg mechanism 
and the dark matter stability is ensured by the choice of quantum numbers. 

• One can generate mass for the new gauge boson using the Stueckelberg or the Higgs mech¬ 
anism. In the Stueckelberg scenario the dark matter candidate interacts only with the new 
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neutral gauge boson in the theory. However, if one uses the Higgs mechanism the dark matter 
could also have interactions with the new physical Higgs boson. 

• The existence of a new gauge boson is key for the testability of the mechanism for dark matter 
stability. Therefore, in the ideal case the new gauge boson should define all the properties 
of the dark matter candidate. In particular, the main annihilation channels must proceed 
through the interaction between the dark matter and the new gauge boson. One can show 
that in this case the dark matter annihilation through the new gauge boson is not velocity 
suppressed. 

For a Dirac dark matter, we investigate in detail the relic density constraints, the predictions 
for direct detection and the dark matter annihilation channels producing monochromatic photons, 
XX —> Z' —>■ 77 , hy, Zy. We compute the one-loop generated vertices Z' 77 , Z'hy and Z'Zy 
needed for the annihilation cross sections. We discuss all the technical details for the computation 
of the loop graphs and we stress the need to check the Ward and Slavnov-Taylor identities to make 
sure the final results are correct. We point out that the effective coupling Z '77 is possible only in 
models where the charged fermions inside the loop have an axial coupling to the Z'. The effective 
couplings Z'h'y and Z' Zy are always present when the Standard Model fermions are charged under 
the new gauge symmetry. Our results can be used for the study of the gamma-ray lines in any 
theory with Dirac dark matter charged under a new gauge symmetry. 

In order to illustrate the main results, we discuss the simplest possible self-consistent dark matter 
model which is most relevant for studying the connection between direct and indirect dark matter 
searches, as its annihilation cross section is not velocity suppressed. In this model, dark matter is 
charged under the B — L gauge symmetry and one has only two annihilation channels into photons, 
XX —>■ Zbl —>■ hy,Zy. We discuss the parameter space for direct detection in agreement with the 
relic density and collider constraints, and we show the experimental limits on the indirect searches 
for the ganrma-ray lines and bb. We show that the interplay between the relic density, collider 
searches and indirect dark matter detection experiments sets non-trivial bounds on these simplified 
models. 


II. SIMPLIFIED MODELS 

We discuss models where the dark matter is a Dirac fermion x charged only under a new gauge 
force. In this context we can understand why the dark matter is stable. For simplicity, we consider 
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the case where one has an Abelian force, i.e., a U{\)'. The part of the Lagrangian relevant for our 
discussion is 

C D Fj l + ixL'fD&L + iXRTfDuXR ~ ( M x xr X l + h.c.) + ^M 2 z ,Z'^, ( 1 ) 

where 


II 

1 

(2) 

D^xl = ( d M + ig'n L Z XL , 

(3) 

D^xr = ( d M + ig'n R Z^ X r- 

(4) 


Here we neglect the kinetic mixing between the new Abelian symmetry and U{\)y- For the moment 
we do not discuss the anomaly cancellation and how the masses are generated but will address these 
issues later in a well-motivated model. Notice that in general the Standard Model fermions can be 
charged under the new symmetry such that new fermions are needed for anomaly cancellation. 

In these models the relevant interaction of the dark matter candidate X = XL + XR t° the Z* 
gauge boson is given by 

-ig'x'y 1 * (n L P L + n R P R ) xZl, ( 5 ) 

where we use the standard projection operators, Pj, = |(1 — y 5 ) and P R = 7(1 + y 5 ). The 
interactions of all other fermions / in the theory, the SM fermions or new fermions needed for 
anomaly cancellation, to the Z' can be parametrized as 

- ig'f 7^ (sy - 9a 7 5 ) /■ z 'r ( 6 ) 

As usual, all charged fermions couple to the photon A M according to their electric charge Qf , 

- (7) 

and the coupling of the fermions to the Standard Model Z can be parametrized as 

where 52 is the SU (2) gauge coupling and 6\y is the Weinberg angle. 

A. B — L Dirac Dark Matter 

The local B — L symmetry is anomaly free once we add three copies of right-handed neutrinos 
to the Standard Model particle content. It is well known that this symmetry could play a major 
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role in neutrino physics. Here we focus on a very simple model with Dirac dark matter charged 
under B — L. The relevant part of the Lagrangian is given by 

£bl A ixf^nX ~ M x xx + \ M z bl Z blA, (9) 

where D^x = + WBL n Z RL ) X, n / ±1, and y = XL + XR- One can generate the gauge boson 

mass through the Higgs mechanism or the Stueckelberg mechanism. Let us discuss both cases here: 

• Stueckelberg Mechanism: The mass of the B — L gauge boson can be generated through the 
Stueckelberg mechanism as discussed in Ref. [7]: 

Cbl A \ (M Zbl Z BLii + d fl a) (M Zbl Z% l + cKV) - (Y v I l Hv r + h.c.) , (10) 

where the gauge transformations are given by 

dZ BL = d^X and 5a = —M Zbl A. (11) 

In this case the neutrinos are Dirac fermions because the B — L symmetry is never broken, 
and the dark matter stability is a result of the choice of the quantum number for the dark 
matter candidate. 

• Higgs Mechanism: One can generate the mass for the B — L gauge boson through the Higgs 
mechanism and at the same time we can generate masses for the SM neutrinos through the 
see-saw mechanism [8-12] via the following interactions: 

- C, v = Y v I l Hv r + \rv r v r Sbl + h.c. (12) 

Here Sbl is a Standard Model singlet and has B — L charge two. Notice that if the B — L 
charge of the new Higgs is different from two, the neutrinos will be Dirac fermions. After the 
U{1)b~l is broken, there is a remnant Z% symmetry which is the reason for the dark matter 
stability. 

These simple models have only four relevant parameters for the dark matter study: the gauge 
coupling the dark matter mass M x , the gauge boson mass M Zbl , and n the B — L charge of the 
dark matter candidate. The relevant interactions needed to compute the dark matter annihilation 
channels are 

-igBLn xi^xZ rl and - WBLn f BL (13) 

where n^ BL is the B — L charge of the Standard Model fermion /. 
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1. Relic Density 


The B — L dark matter candidate x can annihilate into all the Standard Model particles and 
the B — L gauge boson Zbl- Therefore, one can have the annihilation channels 

xx -t qq, vv, ZblZbl 

in both the Stueckelberg and the Higgs scenario discussed above. In the Higgs scenario, one has 
the additional annihilation to right-handed neutrinos. There are two main regimes for our study: 

• M x < Mz bl '■ when the dark matter candidate is lighter than the B — L gauge boson, we 
have the following channels, 

xx ->• z*bl ->• qqM,vv,{u R v R ). 

• Mz bl < M x : when the dark matter candidate is heavier than the B — L gauge boson, one 
has a new open channel which is not velocity suppressed, 

XX ZblZbl- 

In order to test this model at the collider, the invisible decay Zbl —> XX is crucial to establish the 
connection between the existence of the new gauge boson and the dark matter candidate. Therefore, 
we focus on the first regime. 



(14) 


Here n[ is the color factor of the fermion / with mass Mf, s is the square of the center-of-mass 


energy, and T z BL is the total decay width of the Zbl gauge boson. In order to compute the relic 
density we use the analytic approximation [13] 



(15) 


where Mp\ = 1.22 x 10 19 GeV is the Planck scale, g * is the total number of effective relativistic 
degrees of freedom at the time of freeze-out, and the function J(xf ) reads as 



( 16 ) 
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n = 1/3 n = 3 



Figure 1. Dark matter relic density flDM^ 2 vs. the dark matter mass M x for different choices of the dark 
matter B — L quantum number n. In the left panel we use n = 1/3, while in the right panel we use n = 3, 
and we give the relic density for two choices of the mass of the Zbl and the gauge coupling gsL that fulfill 
Mz B l / 9bl = 7TeV. The thin blue band corresponds to the currently allowed relic dark matter density 
measured by the Planck collaboration, Udm^ 2 = 0.1199 ± 0.0027 [14]. 

The thermally averaged annihilation cross section times velocity ( av ) is a function of x = M x /T, 
and is given by 

(av){x) = 5 ^ 2 . t [ ax(s- 4M 2 ) R\ ds, (17) 

SM^K^yx) Jam* \ M x J 

where K±(x) and K 2 (x) are the modified Bessel functions. The freeze-out parameter Xf can be 
computed using 


x f = In 


0.038 g Mpi M x (av)(x 

y/9*Xf 



(18) 


where g is the number of degrees of freedom of the dark matter particle. In Fig. 1 we show the 
numerical predictions for the relic density vs. the dark matter mass for two values of n. In the left 
(right) panel we show the results for n = 1/3(3) when Mz BL /gBL = 7TeV, which is in agreement 
with the collider bounds [15]. As expected, for small values of the gauge coupling one needs to rely 
on the resonance to achieve the right relic density. However, generically one can be far from the 
resonance and in agreement with relic density constraints. 
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2. Direct Detection 


The direct detection constraints must be considered in order to understand which are the allowed 
values of the input parameters in this theory. The elastic spin-independent nucleon-dark matter 
cross section is given by 


<7 


SI _ 

X N ~ 


M nM x ffg L o 

k(M n + M x )*M* bl ’ 


(19) 


where Mjv is the nucleon mass. Notice that a^ N is independent of the matrix elements. The cross 
section can be rewritten as 


<7^v(cm 2 ) = 12.4 x 10 41 (j 




GeV 


lTeV V 

rBL ) 


n 2 cm 2 , 


( 20 ) 


where /i = Mi X M x / (Mjv + M x ) is the reduced mass and vbl = Mz BL /gBL- 

In our case M x Mjy, and using the collider lower bound Mz BL /gBL > 6TeV [15] one finds 
an upper bound on the elastic spin-independent nucleon-dark matter cross section given by 


cr®) y < 9.57 x 10 44 n 2 cm 2 , 


( 21 ) 


for a given value of n. There is also a simple way to find a lower bound on the spin-independent cross 
section. The minimal value of the gauge coupling gsL in agreement with relic density constraints 
corresponds to the case when one sits on the resonance, i.e., Mz BL = 2 M x . Therefore, the lower 
bound on the cross section for a given value of the dark matter mass is given by 

c X N > 7-75 x 10" 42 (3bl) 4 ™ 2 cm2 - (22) 

In Fig. 2 we show the numerical predictions for the direct detection cross section a^ N vs. the 
dark matter mass M x compatible with the relic density constraints. The colored dashed lines show 
the values of a^ N for different choices of Mz BL /gBL compatible with current collider limits. The 
black dash-dotted line shows the minimal direct detection cross section. We show the bounds from 
the LUX [16] and XENONIOO experiments [17], as well as the prospects for XENON1T [18]. One 
can see that for n = 1/3, the scenario for Mz BL /gBL = 6TeV is allowed by the LUX experiment 
for a large part of the parameter space. However, for the case n = 3, the ratio Mz BL /gBL needs 
to be larger than 20TeV in order to satisfy the experimental bounds. Therefore, only the scenario 
when n = 1/3 could be tested at the LHC. Unfortunately, the minimal value of the cross section is 
below the neutrino background [19] and it is very difficult to test this part of the parameter space 
in the current direct detection experiments. 
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n DM h 2 = 0.1199 ± 0.0027, n = 1/3 finM* 2 = 0.1199 ± 0.0027, n = 3 




Figure 2. Direct detection cross section a^ N vs. the dark matter mass M x compatible with the relic density 
constraints for n = 1/3 (left panel) and n = 3 (right panel). The colored dashed lines show for 
different choices of Mz BL /gBL compatible with current collider limits. The black dash-dotted line shows 
the minimal direct detection cross section. We show the LUX [16] and XENON100 [17] constraints, as well 
as the prospects for XENON1T [18]. The orange dash-dotted line shows the coherent neutrino scattering 
background [19]. 


B. Upper Bound on the Dark Matter Mass 

In this section we show that in these simple models it is possible to derive an upper bound on the 
dark matter mass. We focus on the case when the Z' is heavier than the dark matter because only 
then one can test the main properties of these models. The argument is based on the observational 
requirement that the relic density of dark matter produced in the freeze-out must not overdose the 
Universe. Today, we know that Udm h 2 < 0.12 and since the relic density scales as Ddm^ 2 ck 1/(c tv), 
one has a lower bound on the annihilation cross section 

(av) > {crv)o « 3 x 10 -26 cm 3 /s. (23) 

Here (crv)o is the minimal value of the cross section compatible with observations. In order to 
guarantee the validity of a given theory we have to make sure that the maximal value of the cross 
section in the theory obeys the condition (crv) max > ( ctv)q . This is a necessary condition, since 
if it is not fulfilled there is no parameter choice in the model which can make it compatible with 
observations. 































11 


The corresponding cross section has the following structure, 


av(g',M Z ',M x ) = 


ci (g') 4 M 2 


with Y z > = C 2 (g') 2 Mz' ■ 


(24) 


(4 Ml - Mf,) 2 + 

It is obvious that the maximum of this expression given a fixed value of g' is realized when 2 M x = 
Mz'- Examining the functional dependence of the obtained expression one finds that 




Cl 


16 c|M 2 ~ 


3 x 10 _ 2 b cm d /s 


(25) 


In the model with local B — L and n = 1/3, this leads to an upper bound on the mass of 

M x < 12.3 TeV and Mz BL < 24.6 TeV. (26) 

Note that this bound is conservative because at the resonance one would need to perform the full 
average as discussed in Sec. IIA 1. This bound is useful to understand the possibility to test this 
type of model. 


III. GAMMA-RAY LINES 

In this section, we discuss the predictions for gamma-ray lines in detail. First, we give the general 
results for any simplified model with an Abelian gauge symmetry and a corresponding Z 1 , then we 
move on to study numerically the predictions for the minimal B — L model discussed before. 

A. Loop-Induced Couplings 


In order to understand the predictions for the dark matter annihilation into photons, we need 
to compute the loop-induced effective interactions Z 1 77 , Z'hj, and Z'Z'y shown in Figs. 3-5. 

• Z '77 coupling: The coupling (5T//// is generated by a loop of electrically charged fermions 
/ charged also under U{\)', see Fig. 3, and is given by 

p la p2pp\+- s i 


= E 


e^ aa (pi - P 2 )a - Pi a P 2 gp\ + -e V(Taf) Pi aP 2 ^ 


f L 


4 


+ £ [e^p la p2pp u 2 - pi a P2fiP IJ i ) 4 + Y, ^ aP PiaP2piPi +P2Y4 • (27) 

/ / 

Since we are interested in processes with two external photons, one has (p±) = p^Yip'i) = 

£ 

0 and the terms in the second line of Eq. (27) proportional to A^ do not contribute to the 
amplitude. The relevant coefficient functions are given by 
, e 2 Q 2 f g' g^Ni 

4 = - f 4?r 2 [3 + A(«, Mf, Mf) + 2 M 2 0,(0,0, s; M f ,M f , M f )] , (28) 
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Pi 

P - 



Figure 3. Coupling Z '77 generated by a loop of fermions / charged under XJ (l)' and carrying electric charge. 
In the calculation, the crossed diagram has to be taken into account. 


Pi 



Figure 4. Coupling Z'h'y generated by a top loop. In the calculation, the crossed diagram has to be taken 
into account. 


f e 2 Q 2 f g' gf.Nl 

A 7 = f “J* I 2 + A (*. . 


(29) 


Notice that this coupling can be generated only if the axial coupling g A of the fermions in 
the loop to the Z' is different from zero. See Appendix B for the explicit form of the loop 
functions A (s, Mf, Mf) and Co(0,0, s; Mf, Mf, Mf). In the above equations s = (p\ + P 2 ) 2 , 
Mf is the fermion mass, and n[ is the color factor and Qf the electric charge of the fermion. 
We checked that the Ward identities pi^ = STg',^ p 2 v = 0 are satisfied. 


• Z'h'y coupling: In Fig. 4 we show the coupling Z'h'y generated by a top loop. The other 
Standard Model fermions have smaller Yukawa couplings and their contributions are therefore 
negligible. For the top quark g l A = 0, and the coupling is given by 

= C z , hl [2 p^ 2 + (M 2 - s) <T] + p^, hr (30) 

where 

Cz'^ = I 2 * [A (M 2 h ,M„M t ) - A 
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Pi 

P - 



Figure 5. Coupling Z'Z 7 generated by a loop of fermions charged under U(l)' and carrying electric charge. 
In the calculation, the crossed diagram has to be taken into account. 


+ {M 2 - s) [2 + (a + ml - M 2 )C 0 (0, M 2 , a; M U M U M t )\ } . (31) 

Here vq is the vacuum expectation value of the Standard Model Higgs. Notice that for 
processes with an external photon, the Cz’h-y term does not contribute to the amplitude. See 
Appendix B for the explicit form of the loop functions. For this vertex, one can show that 
the Ward identity is satisfied, <5T p 2v = 0. 

Z' Z 7 coupling: The coupling between the Z ', the photon, and the Z can be generated at 
one-loop level as shown in Fig. 5. This coupling can be written as 

g'g 2 eQ f N[ 


cj^fiizcr zf ^ j ^ '< t*[wct 

z ' Zl I67r 2 cos 6\y 


(32) 


where 


= \ e^ a p 2a B{ + B-f \{M 2 - s)e^ a p la + e^p la p 2f) p\ 


+ e M<TQ/3 pi Q P2/3P2^4 + el/aal3 PlaP2l3PiB( 

+ £ uaa ^PlaP 2 /3P 2 Bq + e^ a ^p la p 2/ 3 (piB f 7 +p 2 b{ 


(33) 


For processes with external Z and 7 only, the terms in the second line of Eq. (33) proportional 
to b[ and B~ do not contribute and the relevant B'j functions are given by 


Bf — 
n 2 “ 


+ 


B{ = — 


(M^){( A ^ + B/9 0 

x [3 (M| - a + M|A(M|, M f , M f )) - (2M| + a)A(a, M f , M f )] (34) 

B f g f v {M 2 z - 2 M 2 ) + Afg f A {M 2 + 2 Mf)] {M 2 - s)C 0 { 0, M 2 , a; M f , M f , M f )} , 


(s-M 2 z y 


(A f gH + Bfgfj [sA{M 2 , M f , Mf) - sA (a, M f , M f ) 
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+ (Mf -a) (1 + 2M 2 C 0 (0, M 2 z , s; M f , M f , M f ))] , (35) 

S 6 = 7-^2 + {3M|A(M|, - (2M§ + s)A(s, M f , M f ) 

(s -Ml) v 7 

+ (M 2 - a) [3 + (2 M) + Ml)C 0 (0, M 2 , a; M f , M f , M f )] } , (36) 

B f - _®_ 

7 (* ~ Ml) 

B f 8 = _^ 2 \ 2 ( Af 9 f A + B f g f v) [2M 2 A(Mf, M f , M f ) - (M 2 + a)A(a, M f , M f ) 

[s M z ) 

+{M 2 Z ~ s) (2 + M 2 C 0 (0, Mf, a; M/, M f , M f ))] . (38) 

Here the coefficients A* and B? are given by 

Af = \ (s{ + 9r) and B f = ^ (g f L - g f ^j . (39) 

Using the Ward identity one can write b( as a function of B 3 : 

b{ = -\{s-Ml)B[. (40) 

In the case of a spontaneously broken gauge symmetry, one has to use the Slavnov-Taylor 
identity [20, 21] 


(Afgl + Bfgfy) [A(M 2 , Mf, Mf) - A(a, M f , M f )\ , (37) 


- = M z 6T% Ay , (41) 

which is the generalization of the Ward identity. Here, A is the Goldstone boson. From this 
identity, one finds the relation 

b{ = -M\b{ - i(a - Ml)B{ - 16Bfg f v M 2 f C 0 (0, M 2 , a; Mf, M f , M f ). (42) 

We have checked that these identities are satisfied in our calculations. Since for our case we 
find B~ = 0, one can write 

= - M 2 Z )B{ - 16Bfg f v M 2 f C o (0, M 2 , a; M f ,M f ,M f ) 

= -\{s-Ml)B{-Ci (43) 


It can easily be checked that the following additional relation between the coefficients holds, 

Bf = B f 8 -B[-B f 3 . (44) 

These relations are very useful to cross-check the results and simplify the final expressions 
for the cross sections. We have used Package-X [22] to perform all one-loop calculations and 
have cross-checked the results. 
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Using the above calculations for the loop-induced couplings, we compute the dark matter anni¬ 
hilation cross sections for the different channels: 

• XX ~ 77 : The amplitude for the dark matter annihilation into two photons is given by 


| M(xx -t 77 )| 2 = 


= fe7/ [1 + 2M|C„(0,0,a; M h M h M,)} j . (45) 

The cross section times velocity for this channel in the non-relativistic limit is given by 

2 


£/ K - HI) 

OkR - n L ) 2 ( g’) 2 M 2 s 3 

4Mf, 

(m|, + r|.) 


°{XX -t ll)v = 


£/ ( M l A 7 - l A i) ( n L - n R f (g') 2 M 2 


(46) 


4vr M|, (Mf, +T|,) 

XX h'y: For the dark matter annihilation into the SM Higgs and a photon one finds the 
amplitude 

|2 \ 2 

,~T 


' (»-"!.) + Mf,r|, 


x [cos 2 0 (nf + n|)(s - 4M 2 ) + s{n\ + n 2 R ) + 8 M 2 n L n R \ , (47) 

where 9 is the angle between h and 7 in the center-of-mass system. The corresponding 
annihilation cross section is given by 

(g') 2 \C Z 'h 1 \ 2 (s ~ M 2 ) 3 [( n\ + n 2 R )s - (n 2 L - 6 n L n R + n 2 R )M. 2 ] 


°{xx ^7) = 


487rs 3 / 2 


y/s-uqfa-Aq^ + ri'M*,] 


(48) 


XX Z 7 : In general, the explicit form of the amplitude for the dark matter annihilation 
into a Z and a photon is very involved and cannot be given here. Here we list the result for 
n R = n R = n where the integrated amplitude is given by 

M„ 2 „ 2,„2 

I.A/Uyy —>• Z^v)f z = — 

2ir 


s 


? \M(XX -4 Zjf = < 9 ')71A 


7687T 4 cos 2 9\v 

(2 M 2 + s ) (Mf - S ) 2 (Mf + s) E/ QfNl (Cl + B{M, 


M 2 s [{s-M 2 ,Y + T 2 z ,M 2 ] 

while the cross section is given by 

{.g’Ygl^n 2 ( 2 m 2 + s) (s - Mf ) 3 (Mf + s ) Qf N c (cl + b[m z 
°{XX - 7 - Zi) = — 


(49) 


245767 T 5 cos Mf s 5 / 2 yjs — 4M 2 [(s — Mf ,) 2 + rf,Mf,] 


( 50 ) 






















16 


Notice that these results are very general and can be used for any dark matter model with the 
features discussed above. For simplicity, we show the results for the XX Z 1 only for a vector 
coupling of the Z' to dark matter. In Refs. [23-25], some of these effective couplings have been 
computed and the implications for dark matter models have been investigated in detail. 


B. Gamma-Ray Lines and B — L Symmetry 


In the simple B — L dark matter model one has only vector couplings to the B — L gauge boson 
and there are only two relevant channels for the annihilation into gamma-ray lines: 


XX Z BL h'y, Zj. 


The energy of the line signal for the process XX ~ 7- 7 A" is given by 


E 7 = M x 


1 - 


m x \ 

4 Ml) ’ 


(51) 


where mx is the mass of the particle X. The cross section for the dark matter annihilation into 
the Standard Model Higgs and a photon is given by 

n 2 g 2 BL \Cz BL h~,\ 2 {s-M %) 3 (a + 2M x) 


c(xx -t fry) = 


247rs 3 / 2 


s-4M2 


( s “ M z BL ? + r L, M i 


Zbl Z bl 


(52) 


The explicit expression for the one-loop generated coupling Cz BL h'y can be found using Eq. (31) 
where (jy is replaced by the B — L number 1/3 for the top quark and g' — y gBL- In the non-relativistic 
limit the above cross section times velocity reads as 


2\3 


a (xx ->• hj)v = 


nY BL \Cz BLh ,\ 2 m^-M^ 


32ttM2 


( M Isr “ 4M v) 2 + r Lr M I 


Z B L Z BL 


(53) 


The annihilation cross section for XX z bl I s gi ven by Eq. (50). 

In Fig. 6 we show the allowed parameter space for the indirect detection for the channels XX ^7 
(left panel) and XX z 7 (right panel) compatible with the relic density constraint. We show the 
results in two different scenarios for the dark matter charge n = 3 and n = 1/3. Only in the case 
n = 3 the bounds from Fermi-LAT rule out a small part of the parameter space in the low mass 
region. In both panels we show the full range for the annihilation cross sections (av) xx ^z-y an< I 
( av )xx~^hj which is defined by the collider limits and the resonance regions. The complete region 
between the two curves is allowed in the model. However, the region close to the resonance can be 
ruled out in the near future. 
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n DM h 2 = 0.1199 ± 0.0027 Q D Mh 2 = 0.1199 ± 0.0027 




M x [GeV] M x [GeV] 


Figure 6. Allowed parameter space for the indirect detection for the channels \X b'y (left panel) and 
XX (right panel) compatible with the relic density constraint. The limits of Fermi-LAT [26] and 

H.E.S.S. [27] are shown by the blue and red lines, respectively. See Tab. I in Appendix A for the values of 
the cross section limits. Note that the highest cross section in the indirect detection experiments corresponds 
to the parameters leading to the lowest cross sections in the direct searches. 


In Fig. 7 we show the predictions for the cross section of the dark matter annihilation XX ~A bb. 
Again, we find the full range for the annihilation cross section in agreement with the relic density 
constraints. The bounds from Fermi-LAT rule out a significant fraction of the parameter space in 
both scenarios. As one expects, this bound is much stronger than the bounds from gamma-ray lines. 
It is important to mention that in the resonance region the direct detection limits are irrelevant 
and the best way to probe the model is through indirect detection. 

We have investigated carefully the final state radiation XX f fh which is mediated by the 
Zbl at tree level and contributes to the continuum gamma-ray spectrum. As we have discussed 
above, the ganrma-ray line in this model is possible only at the quantum level and the cross sections 
are much smaller than the final state radiation. Therefore, it is very challenging to distinguish the 
gamma-ray lines from the continuum spectrum. For example, when n = 1/3 and M x = 200 GeV, 
the difference between the ganrma-ray lines and the continuum is only about one percent. Let us 
stress that the possibility to distinguish the gamma-ray lines from the continuum is a requirement 
to actually use experimental limits from line searches to derive bounds on a particular dark matter 
model. In the future one could have a very good energy resolution in experiments such as Gamma- 
400 [28-30] and one can investigate this issue in more detail. 
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n DM h 2 = 0.1199 ±0.0027 



M x [GeV] 

Figure 7. Allowed parameter space for the dark matter annihilation into two bottom quarks compatible 
with the relic density constraint. The experimental bounds from Fermi-LAT [31] are given. 

IV. SUMMARY 

We have studied the relic density constraints and the predictions for the direct detection ex¬ 
periments in simplified models for Dirac dark matter. In these models the dark matter properties 
are defined mainly by the gauge interactions, and one can understand dynamically the dark matter 
stability. We discussed the cases where one uses the Stueckelberg or the Higgs mechanism for gen¬ 
erating the gauge boson mass. We have presented general results for the three possible annihilation 
channels for the Dirac dark matter giving rise to gamma-ray lines, XX “A Hi ^ 7 , -^ 7 - We have 
shown that the channel XX ~A 77 is present only when the Z' has an axial coupling to fermions 
inside the loop generating the effective coupling Z' 77 . The channel XX “A h'y is mainly generated 
by the top quark because today one cannot have heavy chiral fermions which would change the 
Higgs properties. These results can be used for any model with a Dirac dark matter charged under 
a new Abelian force. 

In order to illustrate our results numerically we investigated a simple model based on local B — L 
where the neutrino masses could be generated through the seesaw mechanism and the dark matter 
is a Dirac fermion charged under B—L. In this case there are only two annihilation channels, 
Xx ~A h'y, Z 7 . We have shown the numerical predictions for these channels taking into account the 
relic density and collider constraints. We have investigated the correlation between the Fermi-LAT 
and H.E.S.S. bounds on gamma-ray lines and the annihilation into bottom quarks to show the 
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constraints on these models. The results presented in this paper tell us how much the bounds on 
gamma-ray lines must be improved to be able to rule out or test some well-motivated and simple 
dark matter models. 
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Appendix A: Experimental Constraints on the Annihilation Cross Sections 

Limits on the cross sections are set using the experimental limit on the photon flux, and the 
relations are different for the various annihilation channels considered in this article. 

• XX -1 77 : In the case of the annihilation of a Dirac dark matter into two photons the cross 
section is related to the flux by 

87T 

( av )xx^ll = -^7 ^7 ■ (Al) 

^ann 

• XX ~t X 7 : When one has the annihilation into a particle with mass mx and a photon, the 
relation is 

( <JV )xx~+ x i = f-^'7 \j m x) *^7 > (A2) 

^ann \ v / 

where E x = M x ^1 — has been used. 

With the help of these relations, one can translate the limits from Fermi-LAT [26] and H.E.S.S. [27] 
given for the annihilation into two gammas into the /17 and Z 7 channels. See Tab. I for the 
Fermi-LAT limits in the range 30GeV < < 500 GeV. 

Appendix B: Loop Functions 

In this appendix, we define the loop functions used throughout the article. The function 
A(s, 7711 , 7172 ) contains the logarithmic discontinuity of the Passarino-Veltman Bq function. For 
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Table I. Limits on the gamma-ray lines from the Fermi-LAT collaboration for E 7 > 30 GeV for the R3 region 
with a NFWc dark matter profile [26]. Using the limits on the 77 channel, we derive the corresponding cross 
section limits and dark matter masses for the Z 7 and hj channels. Notice that for the Dirac dark matter 
case the cross section limits have to be multiplied by a factor of two. 


77 channel 

[GeV[ ( av [10 - 29 cm 3 /sl 

h'y channel 

[GeV[ ( av)h -y [10 29 cm 3 /s] 

Z'y channel 

M y [GeV[ (av)z^ [10~ 29 cm 3 /s[ 

31.2 

0.661 

80.4 

8.77 

63.8 

5.53 

33.0 

0.695 

81.5 

8.47 

65.0 

5.39 

34.9 

1.42 

82.7 

15.9 

66.3 

10.2 

36.9 

3.08 

84.0 

31.9 

67.6 

20.7 

39.0 

3.87 

85.3 

37.0 

69.1 

24.3 

41.3 

3.79 

86.8 

33.5 

70.7 

22.2 

43.8 

6.41 

88.5 

52.3 

72.5 

35.1 

46.4 

6.54 

90.2 

49.4 

74.4 

33.6 

49.1 

5.62 

92.0 

39.5 

76.3 

27.2 

52.1 

3.12 

94.1 

20.3 

78.6 

14.2 

55.2 

3.38 

96.2 

20.5 

80.9 

14.5 

58.6 

7.13 

98.6 

40.4 

83.5 

29.0 

62.2 

6.95 

101 

36.8 

86.3 

26.8 

66.0 

4.59 

104 

22.8 

89.3 

16.8 

70.1 

5.18 

107 

24.1 

92.6 

18.1 

74.5 

5.56 

110 

24.4 

96.1 

18.5 

79.2 

3.08 

114 

12.7 

100 

9.82 

84.2 

2.87 

118 

11.2 

104 

8.78 

89.6 

2.87 

122 

10.6 

109 

8.45 

95.4 

2.82 

127 

9.93 

114 

8.01 

102 

5.77 

132 

19.3 

119 

15.8 

108 

5.73 

137 

18.4 

125 

15.3 

115 

15.2 

143 

46.8 

131 

39.4 

123 

15.1 

149 

44.6 

138 

38.0 

131 

10.8 

156 

30.7 

145 

26.6 

140 

5.29 

164 

14.5 

154 

12.7 

150 

10.6 

173 

28.2 

163 

25.0 

160 

8.15 

182 

21.0 

172 

18.9 

171 

13.0 

192 

32.6 

182 

29.6 

183 

6.68 

203 

16.4 

194 

15.0 

196 

13.3 

214 

31.8 

206 

29.4 

210 

9.19 

227 

21.5 

219 

20.1 

225 

13.0 

241 

29.9 

234 

28.1 

241 

18.7 

256 

42.3 

249 

40.0 

259 

10.2 

273 

22.7 

267 

21.6 

276 

41.2 

290 

90.7 

283 

86.8 

294 

41.4 

307 

90.2 

301 

86.7 

321 

17.3 

333 

37.2 

327 

36.0 

345 

15.0 

356 

32.0 

351 

31.0 

367 

48.7 

377 

103 

373 

100 

396 

51.7 

406 

109 

401 

106 

427 

49.6 

436 

103 

432 

101 

462 

39.4 

470 

81.7 

466 

80.3 
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the special case of m\ = m 2 = m it is given by 


A(s, m, m ) 



For special cases, the Passarino-Veltman C$ function can be given as 


Cq( 0, 0, s; m, m, m) = 


and 



Cq(0, M 2 , s; m, m, m) 



(Bl) 


(B2) 


(B3) 
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